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- Abstract - 

O . We study the matter enhancement of T violation in neutrino oscillation with three 

Q . generations. The magnitude of T violation is proportional to Jarlskog factor J. Re- 

cently, the elegant relation, {Am)i2{^m)23{^m)3iJm = A12A23A31J, was derived, 
P^l where Ajj = Am'^-/{2E) and subscript m implies the quantities in matter. Using 

p • this relation, we reconsider how J^ changes as a function of the matter potential a un- 

^ ■ der the approximation |Am^2l ^ I'^^^ial- We show that the number of maxima for J^ 

depends on the magnitude of sin^ 26'i3 and there are two maxima considering the con- 
straint on sin^ 2^13 from the CHOOZ experiment. One maximum of J^ at a = O^A^) 
r> I is given by J/sin26'i2, which leads to the large enhancement of Jm in the case of the 

d • SMA MSW solution. The other maximum at a = 0(Ai3) is IA12/A13I J/sin26'i3, and 

the enhancement is possible, if sin26'i3 is small enough. These maximal values are 
consistent with the results obtained by other methods. 
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1 Introduction 

Solar neutrino experiments have been observing a Ve deficit for a long time [1] 
and the ratio of v^/ve in atmospheric neutrino has implied a v^ deficit [2], which are 
explained by Ve-^^ oscillation and v^-Vr oscillation, respectively. These experiments 
provide strong evidence that there exist masses and mixings in the lepton sector with 
three generations [3]. 

Long baseline experiments [4] and neutrino factories [5] are operated, or planned, 
in order to obtain more convincing evidence for neutrino oscillation. Furthermore it 
could also be possible to observe CP and T violations. 

As the neutrinos pass through the earth in these experiments, matter effects must 
be considered. It has been studied in the context of long baseline experiments [6], and 
in the context of a neutrino factory [7]. T violation is different from CP violation in 
matter and it is pointed out that it is easy to calculate T violation compared with 
CP violation for neutrino oscillation in matter [8]. The T violating part in matter, 
APt- = P(z/q, — > Up) — Pi^up -^ z/q,), (q;,/9 = e,/i, r) is proportional to Jarlskog factor 
Jm [9] of the lepton sector, unlike the CP violating part. The dependence of J^ on the 
matter potential a = y^GpNe is investigated in other works [10, 11]. 

Recently, Harrison and Scott [12] derived the relation 

{^m)l2{^m)23{^m)3lJm = ^12A23A3iJ, (1) 

where Aj^ = A'mjj/{2E) and the quantities with the subscript m are those in matter. 
The inverse of Jm is the square root of a quartic function of a. This means that Jm 
has either one or two local maxima as a function of a. 

In this letter, we present both the exact and approximate form of Jm as a function 
of a using the above relation. It is shown that the number of resonant maxima of Jm 
depends on the magnitude of sin^2^i3. Taking account of the constraint on sin^ 2^13 
from the CHOOZ experiment [13], we show that there exist two maxima. We also 
estimate the maximal values of Jm in the cases of small mixing angle (SMA) and large 
mixing angle (LMA) MSW solutions [14]. 

2 T Violation in Neutrino Oscillation 

We review T violation in three-neutrino oscillations and state the strategy of this 
letter. In vacuum, flavor eigenstates //^(q; = e,fi,T) are related to mass eigenstates 
i/j(i = 1, 2, 3), which have the mass eigenvalues rrii, by the unitary transformation, 

Z/a = UaiUi, (2) 

where Uai is the Maki-Nakagawa-Sakata matrix [3]. The T violating part, AP-rl^'a -^ 

Ufs) = P{ua -^ ^p) — P{.^i3 ~^ ^a), in three generation after traveling a distance L is 

calculated as 

AuL . A23L . A31L 
APt = loJ sm sm sm , [6) 



where 

J ^ \m[U^,U;,U:pp,]. (4) 

In order to obtain the T violating part in matter, we only have to replace Aj^ -^ {Am)ij, 
Tj . ^ (Tj \ . henre I ^ J 

We would like to study the case where large IS.Pt is realized. In eq. (3), APy is a 
product of Jm and trigonometric functions. In the following calculation, we focus on 
the matter effect of Jm which does not depend on L and determine the maxima of Jm- 

As seen in eq. (4), J consists of the product of Uai- It is complicated to calculate 
Jm directly from {Um)ai, which diagonalizes the matter-modified Hamiltonian Hm, 
although the numerical calculation has been performed [10]. However, it is possible to 
calculate Jm without direct calculation of {Um)ai from the relation 

(Am)l2(A„j)23(^m)31<>^m = A12A23A31 J, (5) 

derived by Harrison and Scott [12]. Since the right hand-side of eq. (5) is a constant 
which does not depend on the matter effect, Jm is inversely proportional to a triple 
product of (Am)jj. Therefore, we study the function of the matter potential a such as 

/(a) = [(A„0l2(A™)23(A„)3l]^ (6) 

and determine the minima of /(a). 

3 Triple Product of Mass Square Differences 

In this section, we study the matter effect on /(a). Harrison and Scott [12] suggest 
that /(a) is a quartic function of the matter potential a and in principle its coefficients 
can be written by the parameters Atti?- and Uai in vacuum, although it is complicated 
in practice. We present the exact form of /(a) in relatively simple form by introducing 
new parameters. The coefficients of /(a) is further simplified under the approximation 

|Am^2l ^ I^^^LI- 

First, let us note that (Am)ij included in /(a) are rewritten by the eigenvalues {\m)i 
of the matter-modified Hamiltonian Hm. as {Am)ij = {^m)j ~ {^m)i- The eigenvalues 
{\m)i are the solutions of the equation for t, 

det{Hm -t) = (Ai - t)(A2 - t)(A3 - t) + a{t - 62){t - S,) = 0, (7) 

where 6i{i = 2,3) and Aj(i = 1,2,3) are the eigenvalues of the 2x2 sub matrix Hij{i, j = 
2, 3) and 3x3 matrix H in vacuum. After a calculation, /(a) is expressed as a quartic 
function of a: 

f{a) = [((A„)2-(A„)i)((A„)3-(A„02)((A„)i-(A„)3)]' (8) 

= f4a' + ha' + f2a^ + ha + fo, (9) 

where the coefficients /« (i = 0, . . . , 4) are presented by Aj and Si in the following. 



The coefficients /4 and /o are 

/4 = (^2 - Ss)', (10) 

/o = {(A2-Ai)(A3-A2)(Ai-A3)r. (11) 

By definition (6), /(a) is semi-positive definite, hence, /4, /o > must be satisfied 
taking account of the hmit a ^ oo and a — *> 0. The relations (10) and (11) are 
consistent with these conditions. 
The other coefficients are 

/S = 2[{62 - Ai)(52 - \2){S2 - A3) + (^3 - \l){6s - \2){Ss - A3)] 

-2(62 - b^f\{b2 - Ai) + (^3 - Ai) + (cychc of A,)], (12) 

h = [{S2 - Ai)(53 - A2) + (cychc of X,)f 

-6[i62 - Xi)i52 - \2Ki6s - Ai) + (^3 - \2m6s - A3) + (cychc of A,)], (13) 
/i = 4[(<52 - A2)(53 - A2)(Ai - A3)'(A2 - Ai) + (cychc of A,)] 

+2(A2 - Ai)(A3 - A2)(Ai - X,)[i52 - \i)i6s - A2) + (cychc of A,)], (14) 

which are relatively simple compared with the case where we don't introduce new 
parameters 6i. In section 5, we present the figures using these coefficients. 

Next, let us show that these coefficients are further simplified under the approxi- 
mation IA12I ^ 1^13 1- As 6i are the eigenvalues of submatrix 

where 1 is the unit matrix, they are approximated by 



\TT |2a \TT |2|rr |2a 

e X I K^el ^12 ^ \ I n \TT |2\A i ^^^2 \Ue3\ ^12 /-.^n 

CI2 = Ai + - — rrrr^, d3 = Ai + (1 - \Ue3\ )Ai3 h — — -— — , (16) 



|2\ A I 1*^621 \'^e3\ ^12 



up to the ffist order of A12 using the unitarity condition. Substituting eq. (16) for 6i in 
eqs. (10)~(14) and taking the standard parameterization, f/gi = C12C13, Ue2 = S12C13, 
Ue3 = sisc^*'', where Cij = cos 9ij and Sij = sin%, the coefficients are calculated as 

/4 ^ 4(Ai3)', /3^-24cOs2^l3(Al3)^ /2=^4(Al3)^ 

/l^-24cos2^l2Al2(Al3)^ /o ^ (Al2)'(Al3)^ (17) 

at the leading order. Note that the order of A12 for /j is important when we determine 
the minima of /(a), /i is the ffist order of A12 and /2, /s, /4 are the zeroth order. Its 
difference determines the magnitude of a for each minima. 

4 Matter Enhancement of the Jarlskog Factor 

In this section, we calculate the minima of /(a) using the coefficients (17) in order 
to determine the maxima of Jm- First, we show that the number of minima depends 



on the magnitude of sin^2^i3, and that there are two minima taking account of the 
constraint on sin^ 2^13 from the CHOOZ experiment. Second, we estimate the maximal 
values of Jm and the energies of the neutrino at maxima in the cases of the SMA and 
LMA MSW solutions. 

Let us start with differentiating /(a) in terms of a: 

f{a)' = 4fW + ^ha" + 2ha + f\ = 0. (18) 

Since only /i is 0(Ai2) in fi{i = 1,2,3,4) from eq. (17), in the limit of A12 -^ 0, eq. 
(18) reduces to 

0(4/40^ + 3/30 + 2/2) =0. (19) 

Hence, there exists a solution at a = in this limit. This means that a solution at 
a = 0(Ai2) exists for A12 7^ 0. 

On the other hand, whether another minimum exists or not is determined by the 
discriminant D of the quadratic equation in the parenthesis of eq. (19), 



D = 9fi - 32/4/2 = 4c«3(Ai3)''(l - 9sin^2^ 



13 



(20) 



If sin^ 2^13 > 1/9, there exists only one minimum at a = 0(Ai2) as Fig. 1 (a). If 
sin^2^i3 < 1/9, then there exists another minimum at a = 0(Ai3) as Fig. 1 (b). The 
restriction of the CHOOZ experiment, sin^ 2^13 < 0.10 [13], is included in the case of 
Fig. 1 (b). 




01 0(Ai2) 
(a) Z) < (sin^ 26*13 > 1/9) 




01 0(Ai2) ^^ a 
0(Ai3) 

(b) D > (sin^ 2^13 < 1/9) 



Fig. 1. /(a) has two(one) local minima for D > (D < 0). The CHOOZ experiment favors 
D > and the figure (b). 

(I) The maximal value of J^ at a = 0(Ai2) 
The solution of eq. (18) at a = 0(Ai2) is 



cos 26*12 . 
a = T^^Ai2. 



cos^ 9 



13 



The minimal value is 



/(a) = sin2 2^i2(Ai2)'(A 



13 



(21) 
(22) 



and thus, from eq. (5), the maximum of the ratio is given by 

'J m. r 



J sin 2^12 ' 



(23) 



which is consistent with other works [10, 11]. This means that J^ is largely enhanced 
in the case of the SMA MSW solution. 

We estimate Jm/ J in two MSW solutions: 

.Jrn_{ 12, for SMA MSW, 
J ~\ 1.1, for LMA MSW, ^ ^ 

where we use sin^ 2^12 = 7.2 x IQ-^ (SMA MSW), 0.79 (LMA MSW) [15]. 

The neutrino energy corresponding to the maximum of Jm/ J, from eq. (21), is 

E = ^^^2gi^^^i2 (25) 

where iVg is the electron number density: N^. = 8.2 x lO^^cm^^ in the earth's crust. 
Substituting the experimental data, it is obtained as 

r 25 MeV, for SMA MSW, 

\ 62 MeV, for LMA MSW, ^ ^ 

where we use Am22 = 5.0x10-^ (SMA MSW), 2.7x10"^ (LMA MSW), and sin^ 2^13 = 
0.10(the upper limit of the CHOOZ experiment). 

(II) The maximal value of J^ at a = 0(Ai3) 

The other solutions of eq. (18) at a = 0(Ai3) are 

a = -(3 cos 2^13 ± Vl-9sin2 2^i3) A13, (27) 

where the sign + for cos 2^13 > and the sign — for cos 2^13 < 0. The minimal value 
is 

/(a) = ^^3(^13) [4 _ 3(1 _ 3 sin2 2^13)2 _ cos 2^13(1 - 9 sin' 2^i3)i 1 (28) 

and the maximum of the ratio is given by 



T 



^12 



^13 



1 4^2 , , 

= . (29) 

-13 ' ■ "'- ^ . ■> ^^ -^ ^^ /- „ . .; ^„ X 3 



4-3(1-3 sin^ 2^13)2 - cos 2^13(1 - 9 sin^ 2^13) 



2 



Because of the suppression factor IA12/A13I, the enhancement of Jm is small compared 
with the case (I). 

Furthermore, we can obtain more simple forms for eq. (27) and eq. (29) under the 
approximation 9 sin' 2^13 <^ 1, although this approximation is not justified near the 



upper limit, sin^ 2^13 ~ 0.10, of the CHOOZ experiment. In this case, the value of a 
for the maximum of J™ is 



1 _ sin2 2^13 Ai3 



and the ratio is 



J, 



m 

T 



A 



12 



A 



13 



1 



sin 2^ 



13 



(30) 



(31) 



(32) 



It is understood from this result that the enhancement of Jm is not always realized 
because of the suppression factor IA12/A13I. However, we still have an enhancement 
for small sin 2^13. For example, at sin^ 2^13 = 4.0 x 10^^ which corresponds to sin ^13 = 
1.0 X 10"^ and which is much smaller than the present upper limit, the maximum of 
the ratio is given by 

Jrn_j 0.78, for SMA MSW, 
J ~ \ 4.2, for LMA MSW, 

where we use the experimental values for Am^^ = 3.2 x 10~^eV^. Thus, Jm for the 
LMA MSW solution has an enhancement which is several times as large as J in this 
example. 

The neutrino energy corresponding to the maximum of Jm is calculated from eq. 
(30), 

I sin' 29,,] t"^'' . (33) 



E 



2 ^')2V2GFNe' 

Substituting the same experimental data as before we obtain 

E = 16 GeV. 

We summarize the above two maxima of Jm in Table 1. 



(34) 



a 


Jm/J and E 


SMA MSW 


LMA MSW 


0(Ai2) 


dm/ 'J 


1/ sin 2^12 


12 


1.1 


E 


cos ^Qx^Aml^li^^J^Gpl^ecl-i) 


24MeV 


60MeV 


0(Ai3) 


Jn.lJ 


Am^2/(^^i3 sin 26*13) 


0.78 


4.2 


E 


[1 - (3/2) sin" 2^13] Am23/(2v'2G'FiVe) 


16GeV 


16GeV 



Table 1: The maxima of Jm/J and the neutrino energies E. Input parameters for 
vacuum are the same as in the text, and sin" 26^13 = 4.0 x 10^^ is chosen. 



5 Numerical Estimation of the Ratio Jm/J 

In this section, we numerically study the dependence of the ratio Jm/J on the 
neutrino energy E using eqs. (9)~(15). First, we illustrate the magnitude of maximum 
for Jm taking account of sin" 2^12 and Am"2 given by two MSW solutions and the 
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constraint on sin^2^i3 from the CHOOZ experiment. Second, we study the effect of 
the signs of A12 and/or A13. In the following calculation, input parameters are the 
same as in the previous section and we restrict to the range < % < 7r/4 for simplicity. 
Let us show the energy E dependence of Jm/J in the cases of the SMA and LMA 
MSW solutions with sin ^13 = \Ue2.\ = 0.16 and {U^sl = 1.0 x lO'^ in Fig. 2 (a)~(d) ^ 



(a)SMAMSW:|Ue3|=0.16 



(b)LMAMSW:|Ue3|=0.16 
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(c)SMAMSW:|Ue3|=1.0E-3 



(d) LMA MSW: |Ue3|=1.0E-3 
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1 i—j 1 i_ij 1 i_d 1 . .1 
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Fig. 2. The neutrino energy E dependence of Jm/J in the cases of the SMA and LMA MSW 
solutions with jC/esI = 0.16 and jt/esl = LO x 10^^. The symbol + denotes the maxima 
determined in the previous section. 

Comparing Fig. 2 (a) with (b) (or Fig. 2 (c) with (d)), we conclude that the SMA MSW 
solution has larger enhancement than the LMA MSW solution has for the maximum 
of Jm/J a.t E = O(lOMeV). The enhancement, Jm/J > 1, occurs in the wide energy 
region around this maximum and the values calculated numerically almost coincide 
with the results (24) and (26) obtained approximately. Next, comparing Fig. 2 (a) and 
(c) (or Fig. 2 (b) and (d)), we conclude that if sin ^13 is small enough, the enhancement 
for the maximum of Jm/J Sit E = O(lOGeV) is possible although the energy region is 
small. 

Next, we study the cases where Am^2 and/or Am^g is negative. Since we have 
implicitly assumed that both Am^2 ^ind Am^3 are positive until now, two maxima 
appear in "neutrino" oscillation. However, exactly speaking, whether the maxima 
appear in "neutrino" oscillation or "anti-neutrino" oscillation depends on the signs of 



^Thc energy dependence of Jm/J for a LMA MSW solution with |[/e3 1=0.090 (corresponding to 
Fig. 2(b)) is shown in ref.[10]. 



Aml2 and Amf^. In order to examine such cases, we numerically calculate Jm/J in 
the cases where Amf2 and Am^g respectively, are positive and/or negative, and show 
the results for the SMA MSW solution as an example in Fig. 3. 
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(c) SMA MSW: m3 « ml < m2 
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Fig. 3. The dependence of Jm on the sign of Am?-, (a) and (c) are for Am^2 > , (b) and 



(d) for Ar7T,^2 < 0- i^) ^^^ (b) are for Ayn^g > , (c) and (d) for Am^g < 0. The other 
conditions are the same as in Fig. 2(a). 

Comparing Fig. 3 (a) with (b) (or Fig. 3 (c) with (d)), we conclude that the ap- 
pearance of the maximum for Jm/J dX E = O(lOMeV) depends on the sign of Aml2- 
Although the maximum appears in "neutrino" oscillation in the case of Am^2 > 0, it 
appears in "anti-neutrino" oscillation in the case of Amf2 < 0- Comparing Fig. 3 (a) 
with (c) (or Fig. 3 (b) with (d)), we conclude that the appearance of the maximum 
for Jm/J at -E = O(lOGeV) depends on the sign of Am^g. The maximum appears 
in "neutrino" oscillation in the case of Am^g > 0. On the other hand, it appears in 
"anti- neutrino" oscillation in the case of Am^g < 0. 

These differences originate from the fact that the matter potential a for the max- 
ima of Jm (see eqs. (21) and (27)) is proportional to Amf^-. If Amfj is negative, 
then a is also negative and Jm has the maximum not in "neutrino" oscillation but 
in "anti-neutrino" oscillation. This is because the matter- modified Hamiltonian for 
anti-neutrino is obtained by replacing a -^ —a. 



6 Summary and Discussions 

In this letter, we have studied matter modified Jarlskog factor Jm which appears 
in T violation for the lepton sector in neutrino oscillation. It was shown [12] that the 
inverse of Jm is proportional to the square root of a quartic polynomial of matter po- 
tential a from the relation (Am)i2(Am)23(Am)3i Jm = A12A23A31 J. We have presented 
the exact form of this polynomial with parameters in vacuum and have reconsidered 
the matter enhancement of Jm under the approximation |Amf2l *^ l^^^isl- 

We show that Jm has (i) one maximum at a = 0(Ai2) in the case of sin^ 2^13 > 1/9 
and (ii) two maxima at a = 0(Ai2) and a = 0(Ai3) in the case of sin^ 2^13 < 1/9. 
Considering the constraint on sin^ 2^13 from the CHOOZ experiment, we conclude that 
the case (ii) is realized. 

One maximum of J^ at a = (cos2^i2/cos^ ^i3)Ai2 is given by J/ sin 2^12. Jm/J 
is roughly estimated as 12 for the SMA MSW, thus large enhancement is realized. 
The other maximum at a = (1 — | sin^2^i3)Ai3 is given by IA12/A13I J/ sin2^i3 for 
sin^2^i3 <^ 1/9. If ^13 is small enough, the ratio Jm/J is enhanced. We have roughly 
estimated Jm/J as 4.2 for the LMA MSW solution at sin^ 26*13 = 4.0 x 10"^ 

In the case of sin^ 2^13 = 4.0 x 10^^, our results agree with the results obtained by 
a different method [11]. Our results are also applicable around sin^2^i3 ^ 0.10 which 
is the upper limit from the CHOOZ experiment. 
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